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Abstract 

We investigate the observational consequences of the quintessence field rolling to 
and oscillating near a minimum in its potential, if it happens close to the present 
epoch {z < 0.2). We show that in a class of models, the oscillations lead to a rapid 
growth of the field fluctuations and the gravitational potential on subhorizon scales. 
The growth in the gravitational potential occurs on timescales ^ H^^. This effect is 
present even when the quintessence parameters are chosen to reproduce an expansion 
history consistent with observations. For linearized fluctuations, we flnd that although 
the gravitational potential power spectrum is enhanced in a scale-dependent manner, 
the shape of the dark matter/galaxy power spectrum is not significantly affected. We 
find that the best constraints on such a transition in the quintessence field is provided 
via the integrated Sachs- Wolfe (ISW) effect in the CMB temperature power spectrum. 
Going beyond the linearized regime, the quintessence field can fragment into large, 
localized, long lived excitations (oscillons) with sizes comparable to galaxy clusters; 
this fragmentation could provide additional observational constraints. 

Two quoted "signatures" of modified gravity are a scale-dependent growth of the 
gravitational potential and a difference between the matter power spectrum inferred 
from measurements of lensing and galaxy clustering. Here, both effects are achieved 
by a minimally coupled scalar field in general relativity with a canonical kinetic term. 
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1 Introduction 



Scalar fields are used ubiquitously in cosmology to provide a mechanism for accelerated 
expansion during both the inflationary [H 121 [3] and current dark energy dominated epochs 
(for example, 01 El E])- In the inflationary case, accelerated expansion occurs when the 
inflaton field is rolling slowly and it ends once the field starts oscillating around the minimum 
of its potential. The oscillatory phase is well studied and gives rise to rapid growth of 
inhomogeneities because of couplings to other fields or self interactions (see, e.g. [3 El El HO] ) . 
For a recent review, see [11]. Like the inflaton, the quintessence field provides a mechanism 
for accelerated expansion in the slow-roll regime. However, unlike inflation, the possibility 
of the accelerated expansion ending through quintessence oscillations is rarely discussed in 
the literature (however, see e.g. [12] )• If quintessence was rolling slowly in the past, it is 
"natural" (though certainly not necessary) for it to enter an oscillatory regime as it finds its 
way towards a local minimum in its potential. 

In this paper, we wish to understand the observational signatures of a transition to an 
oscillatory regime in the quintessence field, if such a transition occurs in our recent past {z < 
0.2). An obvious observational signature is a late time change in expansion history. More 
importantly, as is the case with the inflaton, the coherent oscillations of quintessence near 
an anharmonic minimum are unstable and quickly fragment in a spatially inhomogeneous 
manner, giving rise to additional structure. This instability driven by self interactions of the 
field (the anharmonic terms) leads to a growth in structure that is much faster than the usual 
gravitational growth. For example, all potentials that have a quadratic minimum and are 
shallower than quadratic away from the minimum (see figure [T| suffer from this instability 
(in particular for wavenumbers k <^ m, where fn^ ~ U" [ip)\ip^Q. See, for e.g [IS [H El E] ) . 
We find that while the gravitational potential is influenced strongly by such growth in the 
quintessence fluctuations, the over density in baryons and dark matter does not change 
significantly. The scale dependent growth in the gravitational potential and the absence of 
similar growth in the matter power spectrum provides a signature to confirm or rule out 
such a transition in the dynamics of the quintessence field. In addition, if the quintessence 
field becomes nonlinear, it can fragment rapidly into long lived, localized excitations called 

osciiions (see e.g. [ISl [HI [IHl [ISl [201 [211 [221 [2S1 [21 123 12S1 123 12H1 El 12^ Effl with 

size I ~ few x ~ few x Mpc[^ Although we do not pursue this in detail in this paper, in 
appendix B, we comment on how such rapid nonlinear fragmentation of the field, and to a 
lesser extent, the time-dependent evolution of oscillon configurations themselves could have 
interesting ISW signatures. 

While similar resonant phenomenon are commonly invoked in the early universe, for ex- 
ample during preheating {z ^ 10^°) [TT], few observational signatures exist (e.g. [3H1 jSH 
[351 [36]). The reason for this lack of observational signatures on astrophysical scales is that 
resonance is effective only on subhorizon scales. Since high redshifts imply a small Hubble 

*We will frequently relate the mass of the scalar field m, to length and time scales, without explicitly 
writing H and c. We set h = c = 1 throughout the paper. 
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horizon, only small astropliysically inaccessible scales today (of order meters for GUT scale 
inflation) become excited. Moreover, on these small scales, the thermal radiation dominated 
state of the universe at the time of big bang nucleosynthesis smooths out inhomogeneities. 
Since in our scenario the transition in quintessence is happening today, the resonance now 
happens on astrophysically accessible scales. As a result these models can be strongly con- 
strained. 

The observationally interesting phenomenology of quintessence oscillations comes at a 
price. In addition to the usual problems with quintessence models (such as the smallness of 
the energy density, long range forces, etc.), for the quintessence potentials we work with, we 
require that the transition to an oscillatory regime happens close to today. This imposes an 
additional tuning for the initial conditions of the field. 

Related to the present paper [2], Johnson and Kamionkowski discuss the dynamical 
instability arising from a "small" anharmonic term in oscillating dark energy models and 
conclude that such instabilities render oscillating dark energy models unsuitable for pro- 
viding a sustained cosmic acceleration period. While we rely on these very instabilities to 
source the rapid growth of the field fluctuations and the gravitational potential, we envision 
the transition happening late enough so that this instability does not significantly affect the 
observable expansion history. We then calculate and discuss in detail the additional observ- 
ables that can be used to probe the rapid growth in the gravitational potentials arising from 
the dynamical instability. We reiterate, that we are not interested in the oscillatory phase 
providing the accelerated expansion. 

We briefly mention a few papers which explore clustering in quintessence like fields, 
though through very different mechanisms. Observational effects of quintessence clustering 
caused by phenomeno logically varying its speed of sound [not from resonant behavior dis- 
cussed here) has been investigated by several authors (see, for example, [371 EH EHl SOl IHl 
nil Us]). Instabilities in coupled dark energy-dark matter models have also been considered 
in [H]. In addition, rapid transitions in dark energy, in the context of coupled/unified dark 
matter and dark energy models where dark energy "switches on" at late times have been ex- 
plored by |15] (also see references therein). In our case, dark energy "turns off" at late times 
as the quintessence field starts oscillating, leading to a rapid growth in field fluctuations. 
With a somewhat different motivation, in [46j , the authors discuss how rapid, extremely low 
redshift {z < 0.02) transitions in dark energy can be hidden from expansion history mea- 
surements. However, they did not consider the effects of such transitions on perturbations, 
which is the focus of this work. 

In this exploratory paper, we restrict ourselves to single field quintessence models, as- 
sume a minimal coupling to gravity and impose no non-gravitational couplings to other 
fields. We include WIMP (Weakly Interacting Massive Particle) dark matter along with 
the quintessence field and assume a spatially flat universe. We limit ourselves to a linear 
treatment of the fluctuations, including scalar gravitational perturbations, except when we 
discuss the nonlinear fragmentation of the quintessence field and the formation of robust. 
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localized quintessence excitations (oscillons). We focus on a case that is consistent with 
expansion history and galaxy clustering but could potentially be ruled out by measurements 
of the gravitational potential, in particular via the CMB temperature power spectrum at 
large angular scales. Our intention is to stress the interesting phenomenology of such mod- 
els and to point out that such transitions can be better constrained by going beyond the 
measurements of expansion history and growth of structure in galaxies alone. 

The rest of the paper is organized as follows. In section ^ we introduce and motivate 
the form of the quintessence potential used in this paper. In section ([S]) we work through 
the evolution of the quintessence field in an FRW universe. We highlight initial conditions 
and important regimes of evolution as well as constraints placed on the parameters from 
observations of the expansion history. In section Q we discuss the evolution of linearized 
fluctuations in the quintessence field, dark matter and the gravitational potential, with a spe- 
cial emphasis on their evolution during the oscillatory phase of the background quintessence 
field. In the same section we also discuss the domain of validity of our linearized treatment 
and comment on the nonlinear evolution of the field. In section ([s]) we compute observables 
such as lensing power spectra and the ISW contribution to the CMB temperature anisotropy. 
We present our conclusions in section We also include two appendices. In appendix A, 
we provide some details of Floquet analysis and an algorithm used in this paper for calcu- 
lating Floquet exponents. In appendix B, we provide an estimate of the ISW effect resulting 
from evolution of quasi/nonlinear quintessence fluctuations and oscillons. 



2 The model 

We consider a quintessence field governed by a potential of the form (see figure [T]): 



((^/M)2(i- 



where < a < 1. |^ This choice was motivated by monodromy and supergravity models 
of inflation [171 HH HOl |50l [51] and a recent model of axion quintessence [52] (a = 1/2). 
The potential has a quadratic minimum and for very large field values it asymptotes to a 
shallower than quadratic form: 

2 

The scale M determines where the potential changes shape whereas the scale m determines 
the curvature ?7"(0) at the bottom of the well. We will consider cases where the field rolls 



^We write Lp instead of \lp\ to avoid clutter. 
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slow roll 
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Figure 1: Initially, the field rolls slowly causing accelerated expansion of space. As it en- 
ters the oscillatory regime, the acceleration stops and a rapid, exponential growth of field 
fluctuations begins. 



slowly for Lp ^ M, behaving like dark energy and then enters an oscillatory regime with 
(p M after z ~ 0.2. When the field oscillates around the minimum, the "opening up" of 
the potential: 

U{ip) - (l/2)mV^ < ^ 0, (3) 

leads to rapid, scale dependent growth of scalar field fluctuations via parametric resonance 
(see section Q). In the nonlinear regimes, it leads to the formation of localized field excita- 
tions called oscillons [9| [29]. 

Similar phenomenology can arise in models such as pseudo-Nambu-Goldston-Boson (pNGB) 
quintessence (eg. [53j ) 

U{ip) = m'^M^ [1 - cos((^/M)] . (4) 

For potentials ([T]) and (|4]), a limited range of parameters will reproduce the observed expan- 
sion history, allow for a few oscillations in the field close to today and allow for rapid growth 
of structure. For the pNGB like models, the slow roll dynamics of the field necessary for 
accelerated expansion, only occur if M ~ rripi (unless ip — )■ ttM). However, as we discuss 
in section Q, efficient resonance in an expanding universe requires M -C rupi. As a result 
we do not get rapid growth of structure here, except in cases of extreme fine tuning of the 
initial conditions, and do not pursue this model any further in this paper. 

For the potential in ([T]), we find the requirement of a few oscillations close to today 
translates into 

lO^i/o < M,m < IQ-^mp, and a < 1. (5) 
We study these constraints in more detail in the next section. 

Finally, we note that resonant growth of fluctuations also arise when scalar fields oscillate 
in potentials with asymmetric minima 
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and is possible in potentials with a non-quadratic minimum or in potentials which do not 
necessarily open up (though the band of resonant wave-numbers can be narrow). The key 
requirement is that there are anharmonic terms in the potential that provide a time varying, 
periodic frequency in the equation of motion for the fluctuations. 



3 Homogeneous evolution 



In this section, we describe how to choose initial conditions for the background quintessence 
field and parameters in the potential so that the field both oscillates at late times and 
reproduces the observed expansion history. We do this by trying to match the ACDM 
expansion history. This is more restrictive than directly using observations, but makes the 
following discussion more transparent. 

The equations of motion for the homogeneous quintessence field are 



(p + 3i/0 + U'{^) = 0, 



2 



dm 



(6) 



where fldm is the fraction of the critical energy density in WIMPs today, H is the Hubble 
parameter and Hq is its value today The dot "." stands for a derivative with respect to 
cosmic time and the scale factor a = 1 today. We ignore the contribution from radiation 
since we are only interested in late times. 

Current observations are consistent with the quintessence field behaving like a cosmo- 
logical constant in the recent past as well as during the matter dominated epoch [51]. This 
implies that the field's equation of state satisfies w = (0^ — 2U)/{ip^ + 2U) ~ —1 (atleast 
for z > 0.2), which occurs if the 'slow roll' condition 0^ ^ [/ is satisfied. Imposing that 
the energy density in the quintessence field during matter domination has the same order of 
magnitude as the p\ in ACDM, we find: 



(0(S)^°-(-"-'(^)^ 

where (pi is the initial value of the field deep in the matter dominated epoch. Note that for 
simplicity, we've assumed that ^ M and later justify our assumption at the end of the 
analysis. Since we're only interested in scaling relationships, we'll ignore factors of order 
unity from now on. The above gives one constraint between parameters in the potential and 
the field's initial conditions. 

^we choose the subscript "dm" for WIMPs instead of the more usual m to avoid confusion with the mass 
of the scalar field 
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Imposing oscillations in the field at late times gives a second constraint. To derive this, 
we solve equation (|6| in the slow roll, matter dominated regime, assuming the field does not 
change significantly from its initial value. We find (assuming ^ 

^_^^\-a\mlHo)\Mlip,f if « = 0, 

M M ^ y-a^{m/HQf{Mlipif-^"a if a > 0. ^' 

Though the above equation of motion is not valid in the oscillation regime, it gives an 
approximate scale factor (a*) when oscillations begin {^p/M ~ 0). We find: 

{{Ho/mf/\vJM)t if« = 0, 

\(i7o/m)2/3(y,,/M)^«-i ifa>0. ^ ^ 

The actual scale factor when oscillations begin (aosc) differs from a* but the two are mono- 
tonically related. Combining eqs. [7] and |9| we calculate explicit expressions for m and (pi in 
terms of M and a^. We have: 

ifi ja}{mpi/M)^ ifa = 0, 

M \al(mpJM)a-2 if a > 0. 

^ (10) 

m f(mpi/M) ifa = 0, 

^0 ^ \a,^(mp,/M)i-°a-f if a > 0. 

As discussed in more detail below (section |4]), M sets the strength of the resonance 
between the perturbed quintessence field and the background quintessence field. We there- 
fore keep it as a free parameter. We choose so that the actual numerically computed 
oscillations begin around Oosc ~ 0.8. Corrections to distance measurements that depend on 
cosmological parameters beyond Hq are small for scale factors larger than Oosc ~ 0.8. Hence 
modifying the expansion history there will not interfere with measurements that probe the 



cosmological expansion. Last, we fix m and (pi using equation (10). The value of 0, follows 
from equation (|8|. 



^The apparent discontinuity between the cases a = and a > disappears when the fuU expressions are 
used. 
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Figure 2: (a) Homogeneous field evolution: Note that the field is relatively constant until 
ttosc after which it starts oscillating with a decaying amplitude, (b) The Hubble parameter 
increases initially compared to its ACDM counterpart but once the oscillations in the field 
begin, it starts decreasing again. The small oscillations in the Hubble parameter refiect the 
oscillations in the quintessence field, (c) The equation of state parameter w ^ —1 until 
ttosc, refiecting the cosmological constant like behavior due to the slowly rolling field, but 
starts oscillating about w = after Oosc, refiecting the approximately non-relativistic matter 
like behavior of the oscillating field, (d) The acceleration parameter q reveals a period 
of decelerated expansion in the matter dominated era, followed by a period of accelerated 
expansion as the slowly rolling quintessence field begins to dominate the energy density. This 
continues until Oosc, after which q starts oscillating due to the oscillations in the quintessence 
field. 
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In practice, when numerically evolving equation ^ according to the above choices of 
initial conditions and potential parameters, we find that in order to get agreement with the 
ACDM expansion history in the past as well as produce oscillations close to today, we require 
additional fine tuning of the initial conditions and parameters. The procedure described 
above results in a value of Hq that is smaller than the ACDM value. This is because by 
assumption, the Hubble parameter for both ACDM and quintessence are identical. At late 
times, however, the oscillations in the quintessence field cause it to behave like dark matter, 
giving rise to a steeper fall off than the ACDM Hubble parameter. For better agreement 
with Ho, we increase the value of m until Ho agrees to better than a percent. 

Performing the above iteration, we find that a affects whether or not the field starts 
oscillating at low redshift. This should be expected since the asymptotic slope of the potential 
determines the time it takes to transition from slow-roll to oscillations. We can calculate this 
approximate time from the derivative of the quintessence field when slow roll ends. Using 
the slow roll condition described above, the value of the quintessence field when slow roll 
ends is: 

'(mp,/M)i/3M if« = 0, ^^^^ 
arripi if a > 0. 

Assuming 3H(p ~ —U', taking this transition to happen close to today, and using the 



constraints in equation (10), we find the change in scale factor (Aa) between the beginning 



of oscillations {(f ~ 0) and the end of slow roll is given by: 

(aW it« = 0. 

if a > 0. ^ ^ 

The larger the a, the longer it takes for oscillations to start. Hence, consistency with a 
ACDM expansion history and having late time oscillations limits a ^ 1 for the case of 
nonzero a and M <^ rripi for vanishing a. For the rest of this paper, we specialize to a = 0. 
In summary, we keep M and a* as free parameters and choose m, ipi and ipi so that the value 
of Ho agrees with the measured value, the field behaves like dark energy initially, and the 
field starts oscillating at a^sc ~ 0.8. 

For the above prescription, we calculate the evolution of the quintessence field, the Hubble 
parameter normalized by ACDM's Hubble parameter, the field's equation of state parameter 
w and the acceleration parameter q = da/a^ as a function of scale factor. We show the results 
in figure [2] for 

mpi/M = 500, m/Ho = 1130.6, <^i/M = 23.7 with Oosc = 0.82. (13) 

Unless otherwise stated we will use these as fiducial parameter values throughout this paper 
and refer to them as fiducial parameters. As constructed, the field remains approximately 
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constant for a <^ Oosc; rolling slowly towards the minimum. Notice, as assumed initially, 
that (p ^ M aX early times. At Oosc ~ 0.8, the field starts to oscillate with a slowly decaying 
amplitude: 



Note that when Lp ^ M, </9osc('3.) oc a Because of the anharmonic terms in the potential, 
the frequency of oscillation (u) depends on the amplitude of the field and is given by: 



where S is the complete elliptic integral. As — ?■ 0, w — )■ m. Note that in deriving the above 
expression for the period, we've assumed u ^ Hq so that energy is approximately conserved 
during an oscillation. 

The Hubble parameter matches ACDM at early times since both models have the same 
amount of dark matter (see figure [2]^b)). The Hubble parameter must then increase relative 
to ACDM in order for Hq to agree, since the quintessence model's Hubble parameter falls off 
more steeply than ACDM at late times when the field starts to behave like non-relativistic 
matter. The deviation, for the fiducial parameters, is at most 5%. The equation of state 
parameter behaves as expected in the slow roll regime {w = —1) and oscillates between 1 and 
-1 after Oosc as energy swaps between its kinetic and potential parts (figure [2|^c)). Lastly, the 
acceleration parameter reveals a period of decelerated expansion deep in the matter domi- 
nated era, followed by a period of accelerated expansion as the slowly rolling quintessence 
field begins to dominate the energy density. After Oosc, the acceleration parameter starts 
oscillating around q = —1/2 (figure [2|^d)). 

Observational constraints on the expansion history at late times arise primarily from 
distance measurements. Given our above prescription, we expect good agreement between 
distances calculated in our quintessence model and ACDM, which we show explicitly below. 
The comoving distance is given by: 



In figure [3| we plot the percent difference in x(a) between quintessence and ACDM. As 
expected, deviations are smaller close to today, where the Hubble parameter was tuned to 
be identical to the observed Hq, and at early times, where dark energy can be ignored. The 
largest deviations are ~ 4%. Since all observed distances, like the luminosity distance and the 
angular diameter distance, differ from the comoving distance by factors of the scale factor, 
figure [3] implies that current observations cannot distinguish between the background evolu- 
tion of quintessence and the background evolution in ACDM [551 ESI ETJ [581 EH] • We have 



~ V'osc(a) sm{ujt + A). 



(14) 




(15) 




(16) 
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Figure 3: The difference in comoving distance from an observer today, between the 
quintessence and ACDM models. Deviations are smaller today, where both cosmologies 
have the same i^o, and at early times, where both cosmologies have the same amount of 
dark matter and dark energy. 

restricted our analysis to Oosc ~ 0.8 because choosing Oosc ^ 0.8 gives rise to unacceptably 
large deviations from the observed expansion history. 

4 Perturbation evolution 

In this section we investigate the dynamics of the fluctuations in the quintessence field 
{5ip), the gravitational potentials and the overdensity in WIMPs (^dm) assuming linearized 
equations of motion. We ignore radiation and neutrinos, since we are interested in late time 
dynamics. We work in the Newtonian gauge where the metric is 



and $ and \E' are the two scalar gravitational potentials. Tracking the evolution of the grav- 
itational potentials and the two matter components (WIMPs and quintessence) requires two 
second order differential equations. Normally, they are taken to be the energy-momentum 
conservation equations for the matter components. The gravitational potentials can then 
be obtained from the Einstein equations (constraints). However, since the most interesting 
dynamics happen in the gravitational potentials and the field fluctuations, we work with 
them as the degrees of freedom and then obtain the overdensity in WIMPs from the con- 
straints. The equations of motion for the quintessence field fluctuations and the gravitational 
potential are (in Fourier space): 



(1 + 2$)dt2 + ^2(^1 _ 2^)c^a;^ 



(17) 



+ 3/75'(^,. + -^ + U"{^) 5<^fc = -2f/'(^)*fc + 4vj^fc 
■■ • 1 1 r • 1 



(18) 
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The second equation is the diagonal, space-space component of the Einstein equations, where 
the right hand side is the pressure perturbation provided by the scalar field. Since the WIMP 
dark matter is assumed to be pressureless, it does not contribute to this source term. The 
second gravitational potential $a; does not appear in the above equations because there is 
no anisotropic stress at linear order for single, minimally coupled scalar fields. This makes 
the two gravitational potentials equal 



via the off diagonal space-space part of the Einstein equationsj^] After evolving the equations 
for 5ipk and ^dm follows from the time-time components of the Einstein equations: 



-'dm 



01 

ml. 



mi 



mt 



(19) 



4.1 Initial conditions and evolution during matter domination 



We are interested in the behavior of ^v^^, '^^ and 5dm for a > 0.8. To solve (18) we need to 



specify, Sifk-, Sipk, \E'fc, "^k on some initial time slice. This can be done self-consistently using 
the ACDM solutions for \E'fc and 6dm, if the initial timeslice is chosen deep into the matter 
dominated era. During matter domination, the quintessence field is a small fraction of the 
total energy density. As a result, the 6(fk do not contribute significantly to \E'fc and we can 
use from a ACDM cosmology at these early times. We take these initial conditions for 
and directly from the output of CMBfast [SDj at ~ 10~^ since by this scalefactor the 
anisotrotropic stress from neutrinos is negligible and the contribution of dark energy is yet 
to become important. We are then left with specifying the initial conditions in 6(pk, which 
requires understanding its evolution in the matter dominated era. 

The evolution of adiabatic modes on superhorizon: k/aH <^ 1 scales is given by (see 



(H/a) fl a{t')dt' 
-1 + {H/a) Jl a{t')dt' 



The term in the square brackets is constant during matter domination 

'■ 



2^ 
3 H 



(20) 



(21) 



On subhorizon scales {k ^ aH), during the matter dominated era, the gravitational 
potential is determined by the fluctuations in the WIMP overdensity. Since is constant 

Ipor non-minimally coupled scalar fields, — cx f {ip) j f{ip)5ipk where SttGn f{f) ^ ■ In niinimally 
coupled fields, the anisotropic stress is generated by second order terms: $fc — oc 5^\- 
II We will ignore isocurvature modes in this paper 
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during matter domination, the equations of motion ( 18 ) become 



(22) 



0, 



where we have assumed U"{(p)/H'^ ^ {k/aHy. Under the assumption that U'{ip) is slowly 
varying, the complete solution is 



Ck 
1,3 



cos{2kH + 9k) + sin(2A;H + 9k 



2^, 



U'{^) 1 



1 - 



+ 



35 
2k jj 



(23) 



where kn = k/aH and Ck-,9k are constants of integration set by initial conditions at the 
beginning of the matter dominated era. For solutions that are deep inside the horizon, 
kn ^ ^ and only one term grows with time: 



6^k ~ -2^A 



^2 ' 



(24) 



The initial (transient) oscillatory behavior due to the homogeneous part of the solution 23 
as well as the growth oc for -C a ^ aosc due to the particular solution 24 can be 



seen in figure [5|^a). We have also verified that although from equations (21) and (24), we 
see that 6(pk grows as on superhorizon scales and as on small, subhorizon scales, their 
amplitude during matter domination does not become large enough to change the behavior 
of the potential 



With the equations (21) and (24) at hand, we set the initial conditions for all kn using 



a piecewise interpolation between the subhorizon and superhorizon solutions (dashed curve 
in figure [6]). Although we have taken care to faithfully characterize the asymptotic behavior 
of the initial conditions of 6ip, in practice we find that changing the initial conditions by 
orders of magnitude does not affect our results significantly. This is because at late times, 
the particular solution of 6ipk (determined by "^k) dominates, significantly reducing the 
dependence on initial conditions. 



4.2 Evolution during the quintessence dominated era 

In this section, we investigate the dynamics of fluctuations during quintessence domination 
and find that the evolution differs before and after aosc- We analyze these two regimes 
separately. 



4.2.1 Before resonance a < an 



When 0.5 < a < Oosc, the homogeneous quintessence field rolls slowly with an energy density 
larger than, but comparable, to the WIMP density. The behavior of 6(fk on super-horizon 
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scales is still determined by (20), but H decreases more slowly than during the matter 
dominated era. 



In the subhorizon regime, the quintessence perturbations grow faster as we approach a, 



osc 



because the field starts rolling more rapidly (see right hand side of equation (18)). How- 
ever, these fiuctuations are still not strong enough to prevent the decay of the gravitational 
potential, which is caused by a transition to a dark energy dominated epoch. During this 
epoch, the evolution of the '^k is determined by H and does not show any resonant behavior. 
The behavior of the (5^3^ and discussed here can be easily see in figure [Sj^a) and (b) for 
0-5 ^ ^ Oosc- The thin black line represents the evolution of the same mode of in 
ACDM. 



4.2.2 Resonance a > Oosc 

We now come to the most interesting era (a > aosc) with regards to the evolution of fiuctu- 
ations. In this era, the homogeneous field begins to oscillate with an amplitude dependent 
frequency cu < m (see equation ([l5|)). This leads to a rapid growth in the field fiuctuations 
for certain characteristic ranges of wavenumbers. To understand this, we initially ignore 



expansion (H = 0,a = 1) and the gravitational perturbations (\E'fc = 0) in equations (18) 



Resonance in Minkowski space 



The equation of motion for 6(pk then becomes: 

6ipf, + + U"{ip)] 6ipk = 0. (25) 

Since the homogeneous field if is periodic in time, U"{ip) is periodic as well as long as 
U"{lp) 7^ constant, which occurs for potentials with anharmonic terms (as is the case with 
our potential in ([T])). This yields an oscillator with a periodically varying frequency, whose 
solutions can be analyzed via standard Floquet methods (for example, see [62]). For the 
interested reader, we review the main aspects of Floquet analysis in the appendix. Under 
certain conditions (see appendix A), Floquet's theorem guarantees that the general solution 



to equation (25) can be written as 

6Mt) = e^'=*P+(t) + e-^^*P_(t), (26) 



where ±fik are the Floquet exponents and P±(t) are periodic functions with the same period 
as [/"(v^) The Floquet exponent fik depends on the amplitude of the "pump" field yjosc, as 
well as the wavenumber k. There exists an unstable, exponentially growing solution if the 
real part of the Floquet exponent 3f?(/ifc) ^ 0. In figure |4]^a), we show |3fJ(/ifc)| as a function 

**When the Floquet exponents are zero, there exist another class of solutions 5ipk{t) cx tPi{t) and 5ipk{t) oc 
P2{t) where Pi,2{t) are periodic functions. 
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Figure 4: In figure 4(a), we plot the magnitude of the real part of the Floquet exponent 
|3fJ(/ifc)| as a function of the physical wavenumber kp and the amphtude of the background 
field v^osc- Yellow corresponds to a larger value than orange, whereas dark red regions have 
3?(/ifc) = 0. As the universe expands, the wavenumber of a given mode as well as the 
amplitude of the background field redshift. As a result, a given mode "sees" different Floquet 
exponents as it traverses the Floquet chart along the thin white lines. In figure 4(b) we show 
the Floquet exponent "seen" by a mode with a co-moving wavenumber k ~ 0.05m (along 
the left-most dashed white line in Figure 4(a). Modes with k > 0.5m pass through many 
resonance bands in a single oscillation of the homogeneous field and undergo stochastic 
resonance (see text). 

of the amplitude iposc and wavenumber. The color represents the magnitude of the real part 
of the Floquet exponent. Yellow corresponds to a larger |3f?(/ifc)| than orange, whereas red 
corresponds to ^{fik) = 0. Without expansion, neither y^osc nor the momentum k redshift 
with time. As a result, the evolution of modes is determined by the Floquet exponent at 
single point in the (/c,9?osc) plane. This no longer holds true in an expanding universe. 

Resonance in an expanding universe 

The equation of motion for 6ipk in an expanding universe (still ignoring \l/fc) is 



\k^ 



(27) 
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To understand parametric resonance in an expanding background, we make the following 
identifications: ^ 

( ^ -1 (28) 
Here kp is the physical wavenumber and ifosc{0') is the decaying envelope of the oscillating 



field. This identification (28) defines a trajectory in the kp — yjosc plane (white hues in figure 
|4]^a)). The approximate amount of amplification undergone by a given mode is obtained by 
integrating |3fJ(/ifc)| along the corresponding trajectory in the kp — yjosc plane. The |3fJ[/iA;(t)] | 
as seen along such a trajectory [left-most dashed, white line in figure |4|^a)] is shown in figure 
|4|^b). Then, schematically, the evolution envelope of the amplified modes is (ignoring the 
oscillatory piece P±{t)): 



exp 



din a 



H(a) 



(29) 



where in the second equality we use the scalefactor a as a time co-ordinate and 7 = 3/2 
when yjosc <^ M. We assume that the frequency of oscillation u ^ H. This expression 



(29) is meant to give intuition about the resonant behavior in an expanding universe and 
should be used with care, especially when multiple bands are involved since the phase of the 
oscillations can play a role when different bands are traversed. 

When a mode traverses only the first resonance band, the fluctuations get a boost every 
time the homogeneous field crosses zero. This can be thought of as a burst of particle 
production. When a large number of bands are traversed within a single oscillation of the 
homogeneous field, we enter the regime of stochastic resonance. The amplitude of the mode 
still changes dramatically at zero crossings of the homogeneous field, however we are no longer 
guaranteed growth at every such instant. Although over longer time scales the fluctuations 
grow, at some zero crossings they can also decrease. For our scenario stochastic resonance 
is seen for modes with k > 0.5m. A more detailed discussion of these different regimes can 
be found in |63j. 

For our scenario, the modes that grow the fastest are the ones with k <^ m (ones travers- 
ing the first band). To get a sense of what is required of the parameters for these modes to 



grow rapidly, let us now concentrate on the right most expression in (29). If the argument of 
the exponent is significantly larger than one, then we will have rapid growth in fluctuations. 
Using Oosc ~ 0.8, resonance takes place in the logarithmic interval Alna ~ 0.2. Hence, 
one gets rapid growth in fluctuations when |3fJ(/ifc)|/if > 10. During the oscillatory regime 
H ^ Hq m (M/mpi) while from figure |4| the real part of the Floquet exponent seen by 
a mode in the first band has value of |3fJ(/ifc)| ~ 0.1m. Thus we need M/m^i < 10~^ for 
efficient resonance. However, as we saw in section ([s]), in the class of models described by 
([T|, consistency with the observed expansion history and the requirement of a few oscilla- 
tions in the homogeneous field close to today, automatically yields M/rripi < (Aa)^/^ ~ .03 



(equation (12)). Hence in such models, resonance is almost inevitable. 
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Importantly, this is a non-gravitationally driven growth, driven by the homogeneous, 
oscillating field u ^ H. Hence this growth can happen on a time scale which is significantly 
shorter than H~^. 



Resonance in an expanding universe including local gravity 

Let us now include the gravitational potentials in the equations of motion and consider 



resonant phenomenon in the complete coupled system given by equation (18). When the 
field driven resonance is efficient, the gravitational potential does not play a significant 
role in the field dynamics. The evolution of "^k, on the other hand, is affected significantly 
by the resonant growth in S(pk since the quintessence field dominates the energy density. 
We stress that can only be ignored in the evolution of 6ipk for modes undergoing strong 
resonance, when the self interactions dominate over the gravitational one. In particular, as 
we have seen in the previous sections Sipk just before resonance is determined by \E'fc, and as 
a result ignoring before resonance is not justified. In addition, the gravitational potential 
includes contributions from the WIMP dark matter overdensity via the constraint equations. 



These two considerations make the analysis of equation (18) in the resonance regime non- 
trivial and we have to rely on numerical solutions. A more detailed analytical analysis of 
resonance in an expanding universe including the effects of the gravitational potential will 
be pursued elsewhere. Below we discuss the numerical solutions during resonance in a bit 
more detail. 

Typical, rapidly growing solutions for S(fk and with k 0.05m are shown in figure 
and (b) respectively. Note the rapid, nearly exponential growth of \l/fc and S^pk after 
Oosc ~ 0.8 (vertical dashed line). The thin black lines show the evolution of the same modes 
in ACDM. The growth in the gravitational potential is somewhat delayed compared to the 
field. This is because the energy density in the field has to first become comparable to that 
of dark matter. Before this happens, the gravitational potential does not experience scale 
dependent growth (though it still responds to the changing expansion history). 
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Figure 5: (a) Evolution of 5ipk with k ~ 0.05m ^ O.OlMpc^ . During matter domination, 
5ipk a? and grows relatively slowly. However, after aosc, the mode grows exponentially fast 
due to parametric resonance. At a.^\ the fluctuations become nonlinear and linear evolution 
is no longer applicable, (b) Evolution of the scalar gravitational perturbation {k ~ 0.05m). 
Before resonance, the evolution is similar to that expected in usual slow roll quintessence 
(or ACDM, thin black line) model. It is constant deep in the matter dominated era and 
starts decaying as quintessence takes over. Unlike slow roll quintessence, after Oosc, the 
potential grows rapidly until a^i. This leads to a scale dependent signal in observations that 
are sensitive to the gravitational potential, (c) The evolution of the WIMP overdensity is 
not significantly affected by the resonant growth. The departure from ACDM in this case is 
almost entirely due to the slight deviations in the expansion history. The normalization is 
set to one for all the above modes at a = a^. 



So far we have ignored discussing the evolution of the WIMP overdensity since it is deter- 



mined from the constraint equation (19). Naively one might expect to see scale dependent 
departures in the behavior of 5dm- However, despite the rapid growth in the gravitational po- 
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tential after a^sc, ^dm does not deviate significantly from its ACDM counterpart (see|5](c)) in 
the linear regime. Although difficult to see from the constraint equation (19), this behavior 
can be understood by considering the conservation equation for the WIMPs: 



5dm + 2H5a^ = -^^k + 3 (2H^k + ^k) ■ (30) 

Heuristically, we see that on subhorizon scales, 5dm is obtained from a double time integral 
of the potential This delays the response of the (5dm to the \E'fc. The small departure of 
(5dm from its ACDM counterpart can be accounted for by the difference in expansion history 
between the quintessence and the ACDM models. In particular, since the Hubble parameter 
in our quintessence model is always slightly larger than its ACDM counterpart (see section 
[3]), the growth of (5dm in the quintessence model is slightly surpressed. 

From the above discussion we see that the linearized fluctuations in the fleld and the 
gravitational potential grow rapidly. Eventually, the fleld fluctuations will become nonlinear 
and we cannot trust the linearized treatment. Hence it becomes important to understand, 
at least qualitatively, when the fleld becomes nonlinear as well as what happens thereafter. 
Although we discuss the nonlinearity of the scalar fleld fluctuations, we do not include the 
usual nonlinearity in (5dm at late times. This nonlinearity in (5dm is of course well studied, 
but to include it would take us too far beyond the scope of this paper. 



4.2.3 Nonlinearity and oscillon formation 

We cannot ignore the nonlinearity of the fluctuations when U'{ip) ~ U"{ip)5ip. In the oscil- 
latory regime this happen at the scale factor Oni when 

A5<^(fc,anl) ~ V'osc(anl)- (31) 

where A^^(/c, a) = k^Ps^{k, a)/2n'^ with {6{pkSipk') = {211)^ Psip{k)6^ {k — k') is the power spec- 
trum of the fleld fluctuations (see flgure|6](a)). Heuristically, the left hand side characterizes 
the mean square fluctuations over a spatial region of size L, ie. 

= [A,,(A:,a)],^,_, (32) 

In flgure the orange curves show the evolution of the (absolute value of) oscillatory 
homogeneous fleld and its envelope. The black curves show the maximum value of Asip{k, a) 
as a function of a, and its envelope. The location where the two curves intersect Oni 

^ 0.95, 

is taken as the point beyond which the linearized equations cannot be trusted. In flgure [6|^a 



4.1 



we show As^{k,aT^\) (solid orange curve) and As^{k,ai). As mentioned in subsection 
the normalization of Asip{k, ai) is set by A,j,(A;, a^), which we take from WMAP 7 [55] (A^ 
(9/25)A| = 2.42 x 10^^ at A; = 0.002 Mpc"^ and ra, = 0.966). We also note that although 
k ~ 0.05m ^ 0.01 Mpc~^ is where we see the maximum deviation, resonance signiflcantly 
enhances fluctuations for larger k as well, albeit through stochastic resonance. 
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Figure 6: (a) The orange curve shows the rms fluctuations in quintessence {6(p'^) = 
A^(/c,a) = ^Jk^Ps^{k, a)/2Tx'^ at a^i = 0.95, whereas the dashed curve is lO^x the initial 
conditions at = 0.01. Note the large bump due to resonant growth seen in the 6(p field 
at A; ~ 0.05m ^ 0.01 Mpc~^. The higher k regions undergo stochastic resonance, (b) To 
determine when the fluctuations become nonlinear (ani) we compare the rms fluctuations to 
the background fleld. The orange curves denote the homogeneous fleld whereas the black 
curves denote the maximum value of As^{k,a) at each a. Note the enhancement at each 
zero crossing. Once the field fluctuations become nonlinear, different k modes couple and 
the field fragments into localized, long lived energy density configurations (oscillons). 

In order to calculate the effect of the rapid growth of the gravitational potential on 
observables, we need its evolution until today, ie. a ~ 1. However, as discussed above, the 
nonlinearities of the field do not allow us to compute it for a > Oni = 0.95. To remedy this, 
for the linearized calculation, we take a conservative approach and "freeze" the value of the 
potential at ani, setting \E'fc(a) = ^A,.(ani) for a^i < a < 1. More realistically, the potential 
will evolve further in a scale-dependent manner as the scalar field perturbations undergo 
nonlinear evolution. Fragmentation of the scalar field (discussed below) will enhance the 
perturbations, at least on smaller scales (we provide some simple estimates of the ISW effect 
due to this nonlinear evolution in appendix B). Along with the break-down of linearized 
equations, we also note that after Oni, large gradients in the scale field fiuctuations will 
source anisotropic stress, violating our assumption of = $. 

Although in this paper we do not pursue nonlinear evolution of the scalar field fiuctua- 
tions, below we digress slightly and point out some of the interesting phenomenology that 
results. As nonlinearity sets in, the evolution of different k modes becomes coupled, and 
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back-reaction from the perturbations curtails the resonant growth of fluctuations. The ho- 
mogeneous field then fragments rapidly. For the type of potentials considered here (those 
having a quadratic minimum and a shallower than quadratic form away from the mini- 
mum), most of the energy density can eventually end up in localized, oscillatory, long lived 
configurations of the field called oscillons (e.g. [16], [18], |29]). The central density of these 
oscillons can be greater than the background energy density and their sizes are of order a 
few ~ 10~^H^^ ^ 4Mpc. Although oscillons radiate energy through scalar radiation 
(e.g. p!7l l6^ [651 130] ) ■ our limited radial simulations of individual oscillons reveal a lifetime of 
~ 10^m~^ ~ for an oscillon with width of a few and field amplitude of order M. 

This makes them effectively stable compared to current cosmological timescales. Unlike the 
usual oscillons (eg. [121 [IHl ESI [23) which have a relatively stationary energy density profile, 
the oscillons in our model have an energy density that breathes in and out (also see [TJ]) at 
about twice the field oscillation frequency. 

This phenomenon of oscillon production has been studied in the context of the early 
universe, in particular at the end of a similar scalar field driven inflationary period (e.g. 
[321 [ini [HI [I5J). Here, we point out that it will likely happen at the end of the current 
period of cosmic acceleration as well. As with the early universe studies, we expect the 
oscillons to dominate the energy density. However, unlike the early universe, these oscillons 
have sizes which make them astrophysically accessible, with novel signatures such as the 
existence of dark, low redshift clusters with time dependent configurations made of the 
quintessence field. We stress that it is not guaranteed that the universe will become oscillon 
dominated before today. Although for the fiducial parameters, the field fluctuations will 
become nonlinear before today, it can take some time for the nonlinear field to enter a 
state where it is dominated by oscillons. This time can be longer than the time between 
Oni and today. Quantifying this process, requires simulating the full nonlinear dynamics of 
the quintessence field, including nonlinearities of the WIMP overdensity, which is beyond 
the scope of this paper. Given the similarity of the potential with [T3], we expect many 
of the qualitative results to carry over. However, to perform a reliable calculation to be 
compared with observations, one cannot completely ignore the gravitational perturbations, 
which makes the simulations somewhat more challenging. 



5 Power spectra and observables 

With our understanding of the evolution of individual modes, and the validity of the lin- 
earized equations, we now compute the power spectra of the gravitational potential A|,(/c, a) = 
fc^P^r(fc, a)/27r^ and the WIMP overdensity Ps^^{k,a) and their impact on observables such 
as galaxy clustering, lensing and the CMB. 

The calculated power spectra A|,(/c,ani) and Ps^^{k,an\) (with Oni = 0.95) are shown in 
figures [T] (thick orange lines). The thin black lines show the corresponding power spectra for 
the ACDM cosmology. The gray dashed lines show the power spectra at = 10~^. Note 
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Figure 7: (a) Power spectrum of the gravitation potential (orange). The dashed hne denotes 
the initial conditions at = 0.01 whereas the thin black line shows the ACDM power 
spectrum today. Notice the scale dependent growth of the gravitational potential for k ~ 
0.05m = O.OlMpc"^, aside from which the potential decays in a manner similar to its ACDM 
counterpart, (b) The WIMP power spectrum (orange) does not show any features of the scale 
dependent growth, since the WIMPs do not have enough time to respond to the change in 
the gravitational potential. The difference from the ACDM counterpart (black) is primarily 
due to the difference in expansion history between the quintessence and ACDM models. The 
dashed hne represents 10^ x P^.^{k,ai). Note that since we are using 5dm in the Newtonian 
gauge, Pi^^{k,ai) shows an "upward turn" at small k values. 



that since we are using 5dm in the Newtonian gauge, Ps^^{k,ai) shows an "upward turn" at 
small k values. By construction, the initial conditions in the power spectra A*(/c,aj) and 
Pdm(fc, di) in our quintessence model and ACDM agree. The scale dependent departures in 
A,^{k,a) arise after Oosc ~ 0.8. As discussed previously, the enhancement can be trusted 
until Cni = 0.95, after which we freeze the potential until a = 1. The wavenumber where 
the maximum departures are seen is ~ 0.05m ^ 0.01 Mpc~^. For k ^ Hq and k > m, 
the final spectrum closely resembles the ACDM case. For fixed aosc, we have checked that 
the magnitude of departure from ACDM increases with increasing m. This is consistent 



with the notion larger m corresponds to a larger nipi/M (see equation (10)) leading to more 



efficient resonance. As discussed for single modes, the potential A^(k,a) show dramatic, 
scale dependent departures from ACDM but the WIMP overdensity power spectrum does 
not. 

Assuming a model for the bias, the WIMP power spectrum can be probed by measuring 
the two point correlation function of galaxies. Given the lack of scale dependence, and the 
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small difference (< 6%) with ACDM we expect the departures from ACDM to be difficult to 
detect. The fluctuations in quintessence are not directly observable, except through their 
gravitational imprint. This leaves us with the gravitational potential as the key probe to look 
for departures from ACDM. The rapidly changing potential in the quintessence model and the 
large scale dependent departures at low redshifts could affect the low multipoles of the CMB 
temperature and weak lensing power spectrum. With this in mind, we now calculate the weak 
lensing (convergence) power spectra and the CMB temperature anisotropy. We remind the 
reader that the differences that arise here are inspite of having made the expansion history 
consistent with current observations. An expansion history and galaxy power spectrum 
consistent with ACDM but scale dependent deviations in the lensing and CMB power spectra 
provides a unique probe of a late time transition in the quintessence field. 



5.1 Integrated Sachs- Wolfe (ISW) effect 

When CMB photons travel through an evolving gravitational potential, their energy changes 
(the ISW effect [SS])- Since the rapid growth in the quintessence fluctuations induces large 
changes in the potential, we expect imprints of the late time quintessence transition on the 
CMB power spectrum through the ISW effect. 

The angular power spectrum of the CMB temperature can be written as [67] 

Ci = A7r j d\nk Df{k)Al{k)\p,,^ (33) 

where A|,(/c)|prim is the primordial power spectrum. In the above, Di{k) = A/(A;)/\1/^™, 
where Ai{k) on large angular scales (under the assumptions of adiabatic initial conditions 
and instantaneous recombination) is given by: 



Ai{k) = Ar (k) + Aj'''^ (k) 

6 6 CLpilp 



^fc(ae)+ / da ji{kxa)da[^k{a) + <!>k{a)]. 



(34) 



In the above expression, the subscript "e" stands for emission (ie. surface of last scattering), 
Xa = da/a?H{a) is the co-moving distance and ii{x) are spherical Bessel functions. The 
second term is the ISW term, which gets large if the potential evolves significantly between 
last scattering and today. In ACDM the ISW term gets a small contribution just after recom- 
bination since we are not quite matter dominated at that time. But, the main contribution 
on the scales of interest comes from late times as the universe starts becoming A domi- 
nated and the potentials start decaying. The same effect is also present in the quintessence 
model. However, there is an additional contribution from the late time rapid growth of the 
potential. This growth occurs in spurts (see figure [sj^b)) and we can approximate da^k as a 



^^Note that this is reasonable since, hke the WIMPs, galaxies will not have enough time to respond to the 
late time scale dependent quintessence fluctuations. 
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Dirac-delta function at the location of the jumps in the gravitational potential. As a result, 
the contribution of each jump to the ISW term can be easily evaluated: 

(^■)ApW(fc)^2jKA:Xa,)Avl/,(a,), (35) 

where aj is the scale factor where the potential jumps. For a jump of order a few the 
"jump" term is larger than the smooth ISW contribution in ACDM. Note that there can be 
multiple jumps, both positive and negative, depending on the k mode in question. Also note 
that apart from the magnitude of the jump, the location aj also plays a role by fixing the 
argument of the Bessel function. 

We expect that this growth affects the low multipole moments of the CMB since the 
changes in the potential occur very recently. More specifically, the multipole range where we 
expect deviations is / < kxa where k ~ 0.05m ^ O.OlMpc"^. For the fiducial parameters 
chosen here, this yields / < few. Note that along with k, it is the smallness of Xa that limits 
the effect to large angular scales. 

In order to quantitatively calculate the effect of the quintessence perturbations on the 
CMB power spectrum, we use the ^fc(a) and $fc(a) evaluated from CMBFAST [60] for 
a < = 0.01. This captures the contributions from early ISW as well as early anisotropic 
stress (which we ignore at late times). For a > Oj we use \l/fc(a) computed from our own 
independent code for ACDM and the quintessence model. As discussed before, for the 
quintessence case we set da'^k = for a > a^i. With this entire solution at hand (for 
Oe < a < 1) we compute the Ai{k) and Ci. The ratio of the angular CMB power spectra for 
the quintessence and ACDM models is shown in figure |8](a). Note that this is large enough 
to potentially rule out this choice of parameters! Thus, we see that inspite of the effect being 
on large angular scales, ISW provides an excellent probe for constraining the considered 
transition in quintessence. One can change parameters, for example Oosc or m (equivalently 
M) to get the ISW effect small enough so that it is consistent with observations. One 
has to first make sure that such changes are still consistent with the expansion history (see 
discussion in section |3]). Increasing m or aosc or varying m and Oosc in opposite directions can 
lead to an expansion history consistent with observations. Let us now look at perturbations. 
In general, for a fixed Oosc, as we increase m (equivalently decrease M), the fluctuations grow 
more rapidly, shifting the time when the potentials grow rapidly to smaller a values. One 
might expect that this will make the ISW contribution even larger. However, one finds that 
the effect on A/(fc) is not quite as simple. First the change in growth of potential is quite 
sensitive to the number of zero crossing of the homogeneous field. In addition, changing 
parameters changes the wavenumbers which grow the fastest and the time when there is 
a rapid growth in the potential. Heuristically, this affects the argument of the kx of the 
non-monotonic spherical Bessel function in Ai{k). Similar considerations apply to changing 
Oosc- As a result, it is somewhat difficult to apriori predict which combination of parameters 
(consistent with the observed expansion history) will also yield an acceptable ISW term. 
Indeed we find that m ^ 2 x IO^Hq and Oosc ~ 0.82 is entirely consistent with observations 
of the expansion history as well as the CMB. We also find regions of parameter space with 
smaller values of m but larger values of Oosc compared to the fiducial model consistent with 
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Figure 8: (a) Ratio of CMB temperature power spectrum at low multipoles between our 
quintessence and ACDM models. The large difference is due to the rapidly evolving gravi- 
tational potential in the quintessence model after Oosc (b) Ratio of the lensing (convergence) 
power spectrum at low multipoles between our quintessence and ACDM models. The de- 
partures are restricted to low / multipoles due to the proximity of the transition in the 
quintessence field responsible for the rapid growth in the gravitational potential. For lens- 
ing, the model does not agree with ACDM at large / due to the different expansion histories. 

observations. A full sweep of parameters (including a) to determine regions consistent with 
observations is beyond the scope of this paper, but is certainly worth pursuing. 

Before we move on to the calculation of the lensing power spectrum, let us briefly comment 
on the additional ISW that can result from nonlinear field evolution after a^i, ignored so far in 
the linearized treatment. As noted in section 4.2.3, the rapid growth of the linearized scalar 
field perturbations is curtailed once the field perturbations become nonlinear (around Oni). 
However, there will be further nonlinear evolution and rapid fragmentation of the field. To 
get an estimate of the resulting ISW effect, consider a spherical perturbation with mass Ai 
and radius R, collapsing at speed v. Such a perturbation yields a temperature decrement of 
order ~ GAiv/R. For a sphere with an initial radius R^i ~ k~i^ (with k^i ~ 0.05m), an initial 
density comparable to the background energy density and v ~ 0.05 (the approximate group 
velocity of a perturbation with wavenumber kni), the decrement can be ~ few x 10~^ for our 
fiducial set of parameters (m ~ lO^i^o)- The signal is qualitatively similar to the Sunyaev- 
Zel'dovich (SZ) temperature decrement from galaxy clusters [SB]- However unlike < arcmin 
angular scale of the SZ decrement, here, the angular scale is ~ 30 degrees. This simple 
estimate shows that it might be possible to get additional constraints on the quintessence 
transition from the ISW effect due to the the nonlinear field evolution. Individual oscillons. 
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modeled by density configurations with time dependent radii ~ few x can also yield 
an additional ISW signal on smaller scales. We discuss a toy model for estimating the ISW 
signal from such nonlinear perturbations in appendix B. However, we again caution the 
reader that these numbers should be checked with input from detailed simulations of the 
nonlinear field dynamics including gravity. 



5.2 Weak lensing 

Observational constraints from weak lensing observations are often presented in terms of the 
weak lensing convergence (angular) power spectrum [69l [70] : 

= r dxW\x)-^l {k = lx-\x) (36) 
Jo X 



where: 

W{x)= rdx'^^vix), (37) 
Jx X 

A|,(/c,x) is the power spectrum of the gravitational potential and ri{x) is the radial distri- 
bution of sources, normalized to / ri{x)dx = 1- We use the source distribution: 

r]{z) oc ^2exp[-(1.41^/z^ed)'-'] (38) 

with z^ed = 1.26. This distribution approximates the galaxy redshift distribution of the 
COSMOS survey[7T]. As with the ISW effect, because of the closeness of the transition to 
the present day, we expect the signal to be largest at low multipoles / < fcx(ani) ~ few. The 
proximity of the transition also picks out sources at approximately 2x(ani). As a result as 
ttni — )■ l(x — ^ 0), we run out of sources to be lensed. Thus in general, for a fixed enhancement 
of the gravitational potential, we expect the largest lensing signal to arise from the smallest 
Oni still consistent with the expansion history. 

The ratio between the weak lensing convergence power spectrum for the quintessence 
model and ACDM is shown in figure [sj^b). The errors associated with measuring the conver- 
gence power spectrum is given by [69j : 



where /sky is the fraction of the sky covered by the survey, n is the measured galaxy num- 
ber density, and (7int)^'^^ ~ -4 is the galaxy intrinsic rms shear in one component. Using 
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survey parameters characteristic of LSST (/sky ~ .48, n 5.9 x lO^sr"-^) [72], we find 
/\Ci/Ci ~ {.9, .26} for I = {2,30}. Hence this particular quintessence model cannot 
be constrained using LSST measurements of the weak lensing power spectrum. The weak 
lensing convergence power spectrum for the quintessence model does not asymptote to the 
ACDM prediction at higher / because there is a difference in expansion history between the 
two models. Even this deviation would be difficult to see using LSST since the minimum 
ACf /Cf ^ 0.1 at / ~ 300. 

6 Conclusions 

For a slowly rolling quintessence field, it is natural (though not necessary) that the field 
will eventually start oscillating at it approaches a minimum in its potential. In this paper 
we have analyzed the consequences on the expansion history and structure formation of 
a quintessence field, which initially behaves like dark energy and starts oscillating around 
the minimum in its potential at late times (a > 0.8). The potentials considered in detail 
here have a quadratic minimum and are shallower than quadratic away from the minimum. 
When a spatially homogeneous scalar field (if) oscillates about the minimum of an such an 
anharmonic potential, it can pump energy into its spatially inhomogeneous perturbations 
{6(p) through parametric resonance. The amount of energy transfer depends on both the 
wavenumber of 6ip as well as the number of oscillations that take place in the background 
field. This leads to a rapid fragmentation of the homogeneous field and rapid resonant growth 
of the field fiuctuations on time scales significantly shorter than Hq^. 

In order to avoid discrepancies with the measured expansion history, and to simultane- 
ously produce oscillations in the field, we have given a prescription for setting the initial 
conditions of the quintessence field and parameters in the potential. We have also explicitly 
shown that the potentials must be close to constant during the phase where the field is slowly 
rolling. Potentials that don't satisfy this requirement have too much of a delay between the 
end of slow roll and the beginning of oscillations, thus avoiding the rapid resonant growth 
of structure. 

Given a model that produces a background expansion history in good agreement with 
the measured expansion history, we have shown how the gravitational potential and the 
overdensity in WIMPs is affected by the resonant growth of the field fiuctuations. We found 
that the metric perturbations develop scale-dependent growth, with the scale set by the mass 
of the scalar field potential m ~ ^/U"{^ 0) ~ 10^ Ho {k < 0.1m). On the other hand, the 
dark matter overdensity remains featureless and very similar to the ACDM solution aside 
from an overall slight suppression because of the small difference in expansion history between 
ACDM and the quintessence model. This is because, the dark matter doesn't have time to 
respond to the changing potential. Note that scale dependent changes in the gravitational 
potential are normally attributed to modified gravity (e.g. |73l [711 CS]). Here however, we 
have shown that such a change can occur in general relativity with a minimally coupled 
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quintessence field with a canonical kinetic term. 

The rapid growth of the potential significantly affects the ISW contribution to the tem- 
perature angular power spectrum of the CMB. For a range of parameters this yields the 
strongest constraint on such quintessence transitions. Since the metric perturbation devel- 
ops a scale dependent change, the weak lensing power spectrum also offers a possible way 
to constrain this quintessence scenario, where dark energy undergoes a late time transition 
described above. Unfortunately, because of the proximity of the decay in dark energy, de- 
viations from ACDM occur at low / multipole moments. A full treatment, which includes 
nonlinearities in the quintessence field, however, could give rise to deviations in the weak 
lensing and temperature power spectrum, at larger / values. 

The nonlinear dynamics of the quintessence field would give rise to a wealth of new 
phenomenon including field fragmentation and possible formation of localized scalar field 
lumps which could provide additional observational constraints. The full nonlinear analysis, 
combing N-body simulations for the dark matter and lattice simulations for the scalar field, 
is beyond the scope of this paper, but provides a promising avenue to explore quintessence 
transitions and their consequences further. Recall that in spite of its strong clustering 
properties, our quintessence field is minimally coupled with a canonical kinetic term and 
does not have non-gravitational couplings to WIMP dark matter. This could make such 
models easier to simulate than models where gravity is modified (e.g. [761 C3 EH])) the 
quintessence field has non-gravitational couplings to dark matter (e.g. [79]), non-minimally 
coupled quintessence (e.g. |^) or when non-standard kinetic terms are present. 

In summary, if dark energy changes its nature close enough to the present time, it is 
possible to miss it in the expansion history measurements. For the models considered here, 
we have demonstrated that such a transition can dramatically change the gravitational po- 
tential power spectrum in a scale dependent way but leave the galaxy clustering unaffected, 
thus providing a possibly unique signature of such a transition. The best constraints on 
such transitions likely come from the ISW effect, followed by lensing. We expect, a more 
involved analysis, will provide additional constraints on such a transition once we include: 
(i) the nonlinear collapse in WIMPs and in the quintessence field, which will increase power 
on smaller scales (ii) couplings to other fields (ignored here). It would also be interesting to 
explore similar resonant growth in models with multiple ultralight scalar fields (not neces- 
sarily dark energy) motivated in [HI] and studied in more detail in the context of structure 
formation by [52] . 
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8 Appendix A 

Floquet's theorem and calculating Floquet exponents 

The linearized equations of motion for the fluctuations, neglecting the Hubble expansion and 
metric perturbation, are (in Fourier space) 

dld^k + [k^ + U"{^)\ = 0, (40) 

where k is the wavenumber. Since the homogeneous field if is oscillating, U"{(p) is periodic in 
time. This results in a linear system with periodic coefficients, which can be analyzed with 
Floquet theory. Floquet's theorem is most elegantly written in matrix form. Converting our 
second order equation of motion into a first order matrix equation, we find: 

dtx{t) = E{t)x{t), (41) 

where 
and 

Eit) = ^ 1 



~e~u"{^) 

Before stating Floquet's theorem, we need one more definition. The fundamental matrix 



solution, O{t,to), of equation (41) satisfies: 

dtOit,to) = Eit)Oit,to), 

Oito,to) = l. ^ ' 

The fundamental matrix solution evolves the initial conditions x{to) in time: 

x{t) = 0{t,to)x{to). (43) 

Explicitly, O(t,to) consists of 2 columns that represent 2 independent solutions xi,X2, which 
satisfy xi(to) = (1,0) and 0:2(^0) = (0,1)- Since det O(t,to) is the Wronskian of linearly 
independent solutions and det O(to,to) = 1, we have det 0(t, to) = 1. 

We are now ready to state Floquet's theorem (without proof): 
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Consider the linear system 

dtx{t) = E{t)x{t), (44) 

where a; is a column vector and is a real, 2x2 matrix satisfying E{t + T) — E[t) for all 
t. The fundamental solution O{t,to) is 

O{t,to) ^ Pit,to) exp[{t -to)Mito)], (45) 

where P{t+T, to) = P{t, to) and A4(to) satisfies O{to+T, to) = exp[TA^(to)]- The eigenvalues 
/Xi_2 of Ai{to) are called Floquet exponents. 

Since det O{to + T,to) = 1, we have yUi + /i2 = 0. Suppose, from now on, that /ii = 
— /X2 = A*. If M.{to) has two linearly independent eigenvectors e±{to) corresponding to ±/x 
(including /i — 0), then the general solution can be written as 

x{t) = c+V+{t, to)e'^(*-*°) + c_V-{t, to)e-''(*-*°) (46) 

where V±{t,to) — P{t,to)e±{to) are periodic column vectors. If there exists only one eigen- 
vector corresponding to the repeated eigenvalue — 0, then the general solution becomes 

x{t) = c+P+(t, to) + C2 [{t - to)V+{t, to) + Vg{t, to)] (47) 

where Vg{t,to) — P{t,to)eg{to) with eg{to) is a generalized eigenvector. 

From above discussion, we have exponentially growing solutions iff the real part of the 
Floquet exponents ^[fi] ^ 0. In order to calculate Floquet exponents ±fi, it is useful to 
consider the eigenvalues 'jt± of O{to + T, to) called Floquet multipliers. Since det O{to + 
T, to) = 1, then 7r+ • 7r_ = 1. They are related to the Floquet exponents via 7t± = e^'^^. In 
general // and 7r± are complex. Using 7r± = |7r±|e*^='= and 7r+ • 7r_ = 1, we have 

// = ^ [In Itt+I + ie+\ = [In |7r_| + iO^ . 

Thus, 3fJ[/x] ^ iff |7r±| 7^ 1. 

Calculating the Floquet exponents explicitly for the problem at hand reduces to the 
following steps: 

1. First we calculate the period T of U. The period of U (t) depends on the initial amphtude 
of the homogeneous field <^(to) (assuming dt(p{to) — 0) and is given by 

r if max 

n^ma.) = 2 / \ = (48) 

In practice, we specify either (pmax or (pmin- The other is found by solving Vifprnin) — 

y{^max)- For V{ip) = V{-(p) we have ipmax = (Pmvn- 
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2. Next we solve dtO{t, to) = E{t)0{t, to) from to toto + T to obtain: 

O{t,to) = 



6vi\to + T) 6v^i\to + T) 
dtd^^^\to + T) dtS^^p{to + T) 



where {S(p^i}\to 



0}, and {S^f\to) = 0,dt6^'f;>{to) = 1}. This is 



(2), 



equivalent to solving d^cpk + [k"^ + V"{(p)]S(pk — for the above two sets of initial 
conditions from to to to + T. We have suppressed the dependence of T on (fmax to 
reduce clutter. 



3. Last, we find the eigenvalues of O{to + T, to). They are 



(49) 



where all quantities are evaluated at to + T. Note that since the S(pk{to + T) depends 
on the eigenvalues also depend on k. The Floquet exponents ±// are then given by 



±lik = ^ In 



± 



(50) 



9 Appendix B 

ISW effect from nonlinear dynamics of the quintessence field 

In the main body of the text, we discussed how the quintessence field fluctuations grow 
rapidly during the oscillatory regime. This rapid growth leads to a scale dependent growth 
in the gravitational potential which in turn leads to large changes in the CMB temperature 
anisotropics via the integrated Sachs- Wolfe (ISW) effect (see section 5.1). However, after 
Oni ~ 0.95, the field perturbations become nonlinear, and the homogeneous field fragments 
rapidly, potentially forming long-lived, localized excitations of the field called oscillons. Our 
linear analysis did not include the ISW effect resulting from the nonlinear evolution of the 
field perturbations. In this appendix, we remedy this, by estimating the change in the CMB 
temperature due to: 

• The change in the gravitational potential from the initial quasi/nonlinear evolution of 
the field perturbations. 

• The time varying gravitational potential due to an isolated oscillon after it is formed. 
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Let us first look at the ISW effect from the initial evolution of quasi/nonlinear perturba- 
tions. A realistic calculation requires input from detailed simulations of the nonlinear field 
dynamics (for example, [15]). Here we concentrate on a related toy problem of calculating 
the ISW effect due to the collapse of a single spherical top-hat overdensity. This represents a 
crude approximation to the initial rapid quasi/nonlinear evolution of density perturbations 
associated with our model of quintessence. Later, we will evaluate the CMB temperature 
anisotropies from the ISW effect due to a stable oscillon like configuration. Since the universe 
expands very little between ^ 0.95 and a = 1 we will work in a Minkowski background. 
Our main purpose is to obtain an order of magnitude estimate and understand how the effect 
dependence on aspects of nonlinear evolution. 

Consider a spherical mass Ai with uniform density and a time dependent radius R{t) 
located at x^. The gravitational potential due to this mass is given by 



_ GM 
I 1 

*(t,x) 



|x — Xc| > R{t), 



|x-Xc| 

^(^-3) |x-x,|<i?(t). 
The partial time derivative of this potential is 

5t*(i,x) = ^^we(w), (51) 

where u = R^{t) — |x — Xcp and G is the Heavyside function. The change in the CMB 
temperature due to a time varying potential (ISW term) is 

' =2 rft 9t^(t,x)|,=(i„_i)„, (52) 



T 



nl 



where n is a unit vector that defines observer's line of sight, to is the the time since last 
scattering 0.99i/(7^ and we replace x by (to — t)n since light travels on null geodesies. 
We have ignored anisotropic stress here (^P = $). Note that the time of integration starts 
at tni ~ Q.QAHq^ (corresponding to ~ 0.95) since nonlinear perturbations do not exist 
beforehand. 

Let us consider "linear" collapse, ie. 

R{t) = R^x-v{t-U{), (53) 

where < f < 1. This yields 

/I 2\t^ ^^U ^ \ ^ tc-v{Rnl + Vtnl)±rc 

u= {l-v'){t+-t){t-t_), t± = ^ , (54) 

1 — f 

where tc = to — Xc ■ n = to — \^c\ cos 9, Rc = R(tc), Tc = \/ R"^ — {1 — t'^)|xcP sin^ 9 and 9 
is the angle between n and Xc. Note that t_ and t+ denote the time a light ray enters and 
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exist the spherical mass. With these definitions at hand, the integral over time in equation 

&nd 

(55) 



(52) can be evaluated analytically. In particular, if t^i < t- < t+ < to (and Rc > 0) we find 

AT(n) 4GMv(l - v^ fr^^ 



Recalling that and Rc depend on the angle 6, the above expression shows that we should 
expect a cold-spot in the CMB, with an angular size ~ i?c/|xc|. The anisotropy is nega- 
tive and maximal when 6' = 0, when the light ray passes through the center of the mass 
distribution: 

rri I max J-, • W^/ 

I Kc 

The result depends only on the size of the sphere Rc when the light passes through its center, 
the speed of collapse v and the mass A^. To obtain an order of magnitude estimate, we use 
our fiducial set of parameters. For an initial perturbation of size R^i ~ k~^^ !^ (0.05m)~^ !^ 
0.02iJQ'^ (see section 4.2.2 and 4.2.3), and density of order the critical density located at a 
distance |xc| = O.OSHq^ from us (corresponding ~ (to — tni)/2), undergoing collapse with 
V ~ 0.05 (the approximate group velocity of a perturbation with wavenumber /cni), we get 

~ -4 X 10-. (57) 

which could provide additional constraints on our models. This signal is qualitatively similar 
to the Sunyaev Zel'dovich (SZ) temperature decrement from galaxy clusters [HH]- However 
unlike < arcmin angular scale of the SZ decrement, here, the angular scale is ~ 30 degrees. 
As can be checked, the conditions tni < t- < t+ < to and Rc > are satisfied for the above 
mentioned parameters. We remind the reader that this is merely an estimate based on a 
spherical top-hat collapse. The actual nonlinear field evolution can be more complicated. 
In addition, the amplitude of the effect depends on the parameters in a non-trivial manner. 
Nevertheless, this exercise shows that it is indeed interesting to pursue the nonlinear evolution 
of the scalar field perturbations in detail. 

Let us now consider the ISW effect resulting from the time varying potential associated 
with a single oscillon. Based on our numerical solutions (radial only), oscillons with a central 
field amplitude of order Aipjand width of order few m~^, have a peak energy density of order 
m^M^ and an energy density profile that breathes in and out. This is in contrast with a 
fixed energy density profile often used in the oscillon literature. For an oscillon, crudely 
approximated by the a top-hat density configuration, we take 

i?osc(t) ~ 2m-^ ^ i sin^fmt + 0]^ , (58) 

where for our fiducial model m ~ lO^i^Q. We stress that the actual energy density profile of 
oscillons in our model is more complicated. The energy density profile resembles a gaussian 



ttThis M determines the scale where the scalar field potential changes shape, and should not be confused 
with M use to denote the mass of perturbations. 
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when the central field amplitude is at its maximum, but becomes fiatter when the field 
amplitude passes through zero. As a result, a calculation based on realistic oscillon profiles 
could be somewhat different from the estimates below. 

Unlike the linear collapse case, we did not find a simple analytic expression for the 
ISW contribution from individual oscillons. Evaluating it numerically, we find that for our 
fiducial set of parameters and assuming that the energy density within oscillons p ~ im?M'^ ~ 
SHqItiI^, we get 

^^|,nax ~ ±GMoscm X IQ-^ ~ ±10"^ (59) 

We have assumed that the oscillon is located at |xc| ~ O.OSHq^. The shape of the anisotropy 
pattern is ring-like. The phase of the radial oscillation determines whether the center has 
a temperature decrement or increment. Note that the radius associated with an individual 
oscillon is small compared to the quasilinear linear perturbations considered earlier in this 
appendix (by an order of magnitude). As a result, one expects a smaller effect if one assumes 
that their density is still comparable to the average cosmological density. In addition, since 
the light crossing time is comparable to the oscillatory timescale of of the energy density 
configuration, the ISW term undergoes cancellations, leading to a somewhat smaller ISW 
effect than what would be expected on purely dimensional grounds. Thus we expect the ISW 
contribution from individual oscillons to be smaller (and more localized) than that from the 
rapid collapse of nonlinear perturbations. 

Before we end this appendix, we would like to comment on a caveat regarding the emer- 
gence of oscillons. Our linear analysis in the main body of the paper revealed that per- 
turbations become non-linear at a^i- While the non-linear perturbations will form oscillons 
eventually, it is worth asking whether we can see them as large overdensities with individual 
identities by today. To unambiguously understand the timescales associated with emergence 
of oscillons would require a full lattice simulation which is beyond the scope of this current 
paper. Our analysis here indicates that it is certainly worth exploring this further. While 
here, we have concentrated on isolated inhomogeneities, it would be interesting to look at 
the combined effect of a collection of such inhomogeneities (with a number density of order 
(/Cni/27r)^) [9], which would be closer to the actual scenario. 
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